Pseudogap and Amplitude Fluctuations in High Temperature Superconductors 



Philippe Curty^ and Hans Bcck^ 
^ Scimetrica Research, 3007 Bern, Switzerland 
^ Universite de Neuchdtel, 2000 Neuchatel, Switzerland 
(Dated: February 3, 2008) 

Amplitude fluctuations of the pairing field are responsible together with phase fluctuations for 
the pseudogap phenomena in high temperature superconductors. Here we present the more detailed 
theory of the amplitude and phase fluctuations approach in the framework of a fermionic pairing 
model. New experimental comparisons are presented for the specific heat of the curprate LSCO 
confirming the generality of this phenomenological approach. The strong decrease of amplitude 
fiuctuations near optimal doping induces the illusion of a "quantum critical point", which in fact 
does not exist since the pseudogap energy scale is always different from zero even in the overdoped 
regime. 



One major problem concerning high temperature su- 
perconductors is to establish the correct phase diagram 
temperature versus hole or electron doping. A related 
question is the interpretation of the anomalous behaviour 
above the critical temperature. This regime is called 
pseudogap P, because it contains effects similar to 
superconductivity like a partial suppression of electronic 
density of states. The pseudogap region starts below a 
temperature T* where observable quantities deviate from 
Fermi liquid behaviour and seems to be present well in- 
side the superconducting phase. 

Until now two interpretations seem to emerge: the first 
point of view assumes the existence of quantum critical 
point (QCP) sitting in the middle of the superconducting 
regime This quantum critical point is due to an (un- 
known) order parameter which is competition with su- 
perconductivity. The experimental evidences for a QCP 
are however not convincing until now. Moreover effects 
of quantum fluctuatations are suppressed by supercon- 
ductivity at low temperature. Loram and Tallon |j] have 
analyzed many experimental data like specific heat or 
scanning tunneling experiments and extract the value of 
an external energy scale Eg. In their analysis, they as- 
sume the presence of an external gap and extract its value 
by fitting experiments. 

The second group of approaches tends to analyse the 
consequence of pairing and superconductivity on the 
phase diagram [a, 0, H- The pseudogap phase is then 
due to precursor effects of superconductivity or ampli- 
tude pairing fluctuations for example. With this point of 
view, one does not need the existence of the QCP. How- 
ever, there is still the need to explain the characteristic 
pseudogap energy found Eg by Loram and Tallon Q. 

It is perhaps important to remind that T* is a crossover 
temperature that can be directly seen in experiments 
where observable quantities deviate from Fermi liquid be- 
haviour . Eg in the sense of ref. 0| is a phenomenological 
energy scale extracted by fitting experiments according 
to a simple theory ofgapped electrons. 

In a recent letter 0, the importance of treating both 
amplitude and phase fluctuations in a pairing model has 
been shown. Using this method, it is possible to fit accu- 
rately experiments and to extract all the phase diagram 
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FIG. 1: Schematic phase diagram of cuprates. The pseudo- 
gap region of the copper oxides phase diagram lies between 
the critical temperature Tc and a temperature T* . The tem- 
perature Tcj, has been reported in several experiments 0, 
and shows the onset of phase correlations or superconduct- 
ing fluctuations. The energy scale Eg has also been found in 
many experiments like specific heat Q or transport 0. Vb 
and To are the model parameters. Vb is the superfluid stiffness 
and To the mean-fleld pairing temperature. 

of cuprates except the antiferromagnetic regime. By ex- 
panding observables around the average amplitude, one 
can compute thermodynamic quantities like specific heat 
or magnetic spin susceptibility. Two regions have to be 
distinguished (see Fig. : for a relatively small tempera- 
ture interval <T <T^ the phase of ip is still correlated 
in space over some coherence length ^ (the Kosterlitz- 
Thouless coherence length in 2d). Thus, in this regime, 
observables are governed by correlated phase fluctuations 
described by the XY- model. For < T < T*, phases 
of -0 are essentially uncorrelated (^ is on the order of 
the lattice constant), but is still fluctuating and non- 
zero, signaling local pair fluctuations. This explains the 
wide hump between and temperature T* seen in spe- 
cific heat experiments , the depression of the spin sus- 
ceptibility |l3 and the persistence of the pseudogap for 
T < T*. 

One important result is the extracted phase diagram: 
by fitting the experiments it is possible to get the phase 
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stiffness Vq, the mean- field temperature Tq and the 
energy scale Eg versus doping. 

In this paper, we will show that the phase diagram of 
cuprates has indeed an energy scale Eg simililar to the 
one found in 0| . However this pseudogap energy scale is 
due the presence of a pairing attraction and controlled 
by amplitude fluctuations. When these two effects are 
taken into account, they produced a very large Eg in the 
underdoped regime where the superfluid density is low 
and fluctuations are large. In the maximum Tc regime, 
Eg becomes rapidely small because amplitude fluctua- 
tions decrease rapidely. The interaction strenth is also 
decreasing when overdoping. This creates the illusion of 
a QCP in the middle of the superconducting region. 

Our approach has the advantage that the only strong 
assumption is that high temperature superconductivity 
is caused by a fermionic pairing attraction which seems 
to be likely the case. All other assumptions are only 
technical and related to the calculations. 

Various experimental observations can indeed be in- 
terpreted in terms of fluctuations of the pairing fleld 
'0 = iV'le*'^, and that two temperature regions have to 
be distinguished (see Fig. ??): for a relatively small 
temperature interval Tc < T < the phase of is still 
correlated in space over some coherence length ^ (the 
Kosterlitz-Thouless coherence length in 2d) whereas the 
amplitude \ip\ is almost constant. Thus, in this regime, 
observables are governed by correlated phase fluctuations 
described by the XY- model. For < T < T* , phases 
of ip are essentially uncorrelated is on the order of the 
lattice constant), but |?/'| is still non-zero, signaling inde- 
pendent fluctuating local pairs. This explains the wide 
hump seen in specific heat experiments , the depression 
of the spin susceptibility and the persistence of the 
pseudogap for T < T*. Moreover, a magnetic field that 
destroys this pseudogap has to break fluctuating pairs 
and must therefore be much higher than the one which 
suppresses phase coherence and thus superconductivity 

a. 

Our approach has a major difference with the Emery 
and Kivelson phase fluctuations scenario of the pseu- 
dogap regime: our calculations show that phase fluctua- 
tions influence the pseudogap only up to a temperature 
which is much smaller than T*. Above T^, observ- 
ables arc thus only determined by the amplitude of the 
pairing fleld. 

The existence of a temperature intermediate between 
Tc and T* has also been mentioned by Dcvillard and 
Ranninger p^ : using a Boson- Fermion description of 
pairing in cuprates, they flnd that uncorrelated pair- 
ing of electrons leads to the opening of a pseudogap 
at T* . These pairs acquire well behaved itinerant fea- 
tures at T^, leading to partial Meissner screening, and 
thus to diamagnetic susceptibility, and Drude-type be- 
haviour of the optical conductivity. As a function of lat- 
tice anisotropy (and thus of doping) T^ has the same 
tendency as Tc, whereas the higher temperature T* has 




FIG. 2: Amplitude of the pairing field above Tc on a 80 x 80 
array by cluster Monte Carlo simulations: dark color repre- 
sents large amplitude and light colors low amplitude. 

the opposite trend. Although in ref ^3 is related to 
the temperature where the pair life time becomes long, 
it could be identifled with our T^. 

I. MODEL AND EFFECTIVE ACTION 

We are not primarily interested in the origin of pairing 
between electrons but in the consenquences of pairing 
on the phase diagram or physical quantities like specific 
heat or magnetic susceptibility. We assume that there 
is an attractive force between fermions, and we base our 
calculations on a d-wave attractive Hubbard model 

with a hopping t between nearest neighbour sites i and j 
on a square lattice. The interaction favours the formation 
of onsite c?-wave pairs since 

Q^w^E^'^^l (2) 

j 

where Dij = 1,(— 1) for i being the nearest neighbour 
site of j in horizontal (vertical) direction. 

Qh-{4Ai-4A,)/V2 (3) 

is an operator creating a singlet pair on neighbour- 
ing sites. Decoupling the interaction with the help of 
a Stratonovich-Hubbard transformation, the partition 
function Z = Ti e~^^ is then 

Z = Z,J D\,(Te-iy^^'^^^^^'+^ <3lW+hc)^ 

"""(4) 

where 

^ = ^(z,r) = |^(i,r)|e^^(''^), (5) 
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and Hn = — j)a * '^L'^jo- is the non-interacting part. 
The trace over the fermionic operators can be evaluated 
yielding 



Z 



'^''-E.MV'l^+TrinG] 



(6) 



Here G is a Nambu matrix of one-electron Green func- 
tions for fcrmions interacting with a given, space and 
time dependent pairing field ip{i,T). The Green func- 
tions are solution of Gorkov's equations (see 

Expanding ^ in power of V?/;, Z can be written as 
a functional integral involving an action S[ip] for a field 
tp that changes slowly in space and that can be taken 
time- independent : 

^M = ^o(|^|) + ^i(W) (7) 
where Sq is a local functional of ip: 

SoM) = ^ ^ - |^log[2cosh/3£;,/2] (8) 

and S'l = c/ (i'V| VV'P/2 can be considered as the 
deformation or kinetic energy where c is a constant. The 
quasi-particle energy is Eq = + i^q ^ m)^- 

Now we would like to compute thermodynamic observ- 
ables such as energy U = {S)s, specific heat and spin 
susceptibility. The point is that we want to keep the 
XY universality class of the transition together with the 
fermionic character of the system: in the limit of high 
density or weak interaction, the superconductor should 
be described by a BCS like mean field theory whereas in 
the low density limit with strong interaction the transi- 
tion becomes XY like. Our goal is to derive a theory that 
describes these two regimes and, of course, the interme- 
diate regime. 

Our main strategy will be to neglect amplitude 
correlations since simulations show that they are weak 
between different sites — (IV'P) ~ 

since the amplitude is always positive and cannot show 
any critical behaviour. In this spirit, two different 
approaches are possible: 

First approach: the amplitude is fixed but still 
temperature dependent, and is determined by a suitable 
variational equation. Fluctuations from the XY model 
are kept as well as the amplitude weight coming from 
the Jacobian of the cartesian to polar coordinate trans- 
formation. 



II. VARIATIONAL METHOD 

The integration of the partition function can be ex- 
pressed in polar coordinates using the transformation: 

/+00 r + oo p27T 

I -°° i -^0 "'0 

(9) 

Rewriting the free energy F in terms of a constant am- 
plitude \'4>\ yields 

F = _]_iogf £)<^ e-'^('^«(''^"-'°s(l'^l)^/'3+'5'i) (10) 
P J 

where the Jacobian \ipi\ of the polar transformation is 
put into the exponential, and V is the volume. Taking 
the derivative of F with respect to \tp\ and equating it to 
zero leads to the self-consistent equation: 



aso(l*l) V 



Evaluating equation Hll|l and multiplying it with |i/'|/2 
r^-^^ tanh (v/e + H\y{2T)) 



1 

'2P 



iVe 



i<h\2 



0. (12) 



The first and second terms arc the amplitude contribu- 
tion and leads to the BCS gap equation if other contribu- 
tions are neglected. The third comes from the Jacobian 
and implies that the amplitude is never zero. The last 
term is the expectation value of the energy Uxy in the 
XY model with a constant dimensionless coupling K 



K = 



(13) 



where Vb is the zero temperature phase stiffness, and 
is the zero temperature amplitude . This contribu- 
tion characterises the influence of the phase fluctuations. 
Uxy{K) is a monotonic decreasing function with an in- 
flexion point at Tc. Solutions of equation ((IJ) are rehable 
for all temperatures except for T « Tc- However they 
are only expected to be accurate at Tc if the average 
amplitude is large and not varying to much with temper- 
ature. 



III. AVERAGE VALUE METHOD 



Second approach: the energy is expanded around 
the average amplitude. Higher powers of amplitude 
fluctuations are neglected. Here the local coupling 
between phase and amplitude is kept, and amplitudes 
are allowed to fluctuate. 



Since simulations show that amplitude correlations are 
weak, we can expand the energy around the average am- 
plitude where it is not coupled to the phase. The energy 
is then 

U^{S)s^So{{\m + {Si)s (14) 
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FIG. 3: Dimensionless energy of the XY model as a function 
of the coupling K defined in equation H13^ . The continuous 
line is for 3 dimensions and dashed is for 2 dimensions. Note 
that the phase transition occurs near the cusps. 



FIG. 4; Average amplitude (thick line) and the corresponding 
standard deviation (dashed line). Parameters are Vb = 215, 
To = 140. Note that Tc is about 80K whereas the amplitude 
remains non zero at least up to 300K. 



Here, S'o((|V'|)) is the first term of an expansion of the 
average (S'o(lV'l)) around 



5o((|V|))+ci<5|7^|+C2 



^o((H)) + ES(^I^i'") 



(15) 



where 5\il)\ = — . Assuming that amplitudes have 
gaussian fluctuations around their average value, we can 
use the Wick's theorem: 

m\~{mfn = {'2n+im\^?)-m?T 

This last identity together with equation H15|l leads to 



(2ri+l)!! 



2\n 



(16) 



n=l 



Higher corrections term are proportionnal to powers of 
the square of the standard deviation (1-01^) — (IV'I)^- Al- 
though amplitude fluctuations are large, the standard 
deviation is half of the average amplitude, powers of 
(IV'I^) ~ (IV'I)^ s-re small compared to the average am- 
plitude itself. This is why we only take the first term of 
the expansion which already contains the main thermo- 
dynamic features. Additional terms would just add more 
fiuctuations and small corrections to the results. 

For simplicity, averages are computed u sing a nor- 
mafised Ginzburg-Landau action Sgl (see [ij) whose 
potential part Ugl is equal to the first two terms of the 
expansion of So with respect to /SjV'h 



where 



SGL['>P]^kBV^ / d\{UGL + Si) 



UGL=7f{t\^? 



(17) 



(18) 
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FIG. 5: Amplitude distribution for T = 300K (thick line) 
and T = 30K (dashed line). Parameters are Vo = 215K, 
To — 140K. Note the broadening of the amplitude distribution 
at high temperature, i.e the increase of fiuctuations, compared 
to the low temperature distribution where the distribution is 
narrower. 



and t « T/Tq — 1 is the reduced temperature, ip = 
tl'/\ip{T = 0)1 is the reduced field, Tq is the mean field 
BCS pairing temperature. Sgl is normalised with a lat- 
tice spacing e. 

where is the mean field coherence length at zero tem- 
perature, and Vq is the zero temperature phase stiffness. 
Contrary to the variational method, amplitude and phase 
are still locally coupled through Si . The energy becomes 



u^SoiM) 



GL) 



{S 



1)GL 



(19) 



where 



SoUm) = {V {\i'\)VU) - -5^1og[2cosh/3i?,/2] 

9 
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corresponds the BCS free energy for which the gap value 
is determined by the GL average. The quasi-particle en- 
ergy is Eg = [{sq — /i)^ + (IV'I)^ cos^(20)]^/^ where fi is 
the chemical potential. The d-wave symmetry manifests 
itself by the angle dependent amplitude \tp\ cos(20) where 
9 is the angle in q space with respect to qx direction. 

The value of rj depends on the corse-graining procedure 
and is fixed for each sample. Obscrvablcs arc not very 
sensitive to changes in rj. 

Both approaches are valid below and above the criti- 
cal temperature Tc which is the temperature where the 
phase stiffness becomes zero. However the average value 
method gives good results for all values of Tq and Vb 
whereas the variational method works better in the un- 
der doped regime, i.e. for Vq < Tq. 

It is important to notice that the amplitude is fluctu- 
ating although it is fixed to its averaged value in obscrv- 
ablcs. In figureEl the average amplitude and its standard 
deviation are shown. One can note that the latter is ap- 
proximately half of the averaged amplitude. Since the 
amplitude distribution around its average value is almost 
gaussian, the probability density has the following 

form: 

,, -dv-i-d-^i))^ , , 

pdV-D^e TRTR (20) 

This means that the amplitude has values ranging from 
to 2 (I?/) I) causing the large average value 

Computer simulations of the statistical ensemble {"0} 
under the action Sgl have been done using a standard 
Monte Carlo procedure to update amplitude \ip\ and a 
Wolff ^3 algorithm for the phase (p in the same way 
as for the real model O. Typically 10'' sweeps are 
needed to obtain good statistics. . 

Since the fitting procedure is done completely auto- 
matically, one has to simulate first the Ginzburg-Landau 
action ((T7|l for all possible parameters Vq and Tq. Once 
a mesh of simulations has been done, it is possible to in- 
terpolate the surface defined by (|?/'|)(Vo, Tq) and to have 
access all values of (|'0|)(Vb, Tq). The same is done for the 
energy (|VV'P)(Vb, To)- This allows to have an efficient, 
automatic and reproducible fitting procedure. 

IV. SPECIFIC HEAT 

In both approaches, the specific heat C is the sum 
of the amplitude Co and the gradient Ci contributions. 
Defining the reduced specific heat 7 = C/(7„r), we have 

7 = 70 + 71 (21) 

where 71 is divided by Tc instead of T since 5*1 is classical 
and does not satisfy the third law of thermodynamics. 

In practice, the amplitude specific heat can be calcu- 
lated by using the entropy. Co is proportional to the 
derivative of the entropy times the temperature: 

Co[T,\^p\) ^ -T — (22) 




50re 100 150 To 250 r 
Temperature [K] 



FIG. 6: The reduced specific heat 7 from the variational 
method (thick) , which is the sum of the gradient 71 (dashed) 
and amplitude 70 (dotted-dashed) contributions, reproduces 
measurements of YBa2Cu306. 73 (points). Inset: The dashed 
line is the temperature dependent amplitude from equa- 
tion iHTi . 

where the amplitude {ipl is replaced by {ipl in the vari- 
ational approach and by {\tp\) in the average value ap- 
proach. The entropy S for a fermionic system in the 
presence of a gap is a universal function of the ratio 
H'\/T: 

5(|^|/T) = -fcs^/Uog(A.) (23) 

k 

where the Fermi distribution is ft. = — -rn^ — ; — and the 

energy Eh is Ek = £,k + IV'P- By using the entropy, it 
is not necessary to perform the sum over k each time one 
wants to evaluate the specific heat. It is then sufficient 
to take the derivative of the entropy. 

Normalisation of the specific heat: 70 is 1 at high 
temperature since it is divided by the Sommerfeld con- 
stant 7„. The phase contribution is normalised as: 

Ci ^ _k^C^ ^24) 

7n ?o7n 

where Ci is the specfic heat per volume V = N^^q, and 
C^^ /{NkB) is the specific heat per number of lattice sites 
coming from the simulations. Experiments give 7„ ~ 
26 mJ mol-i = 252 J K"! m'^. For the fit of 

Fig. El using the reasonable value ^0 ~ 16 A, we get the 
dimcnsionlcss constant a = ks / [Coin) ~ 13.5. 

A. Specific heat for the variational method 
(d-wave) 

Using the variational method the specific heat C is now 
7 = 7o(K"|)+7i (25) 

where the amplitude \ip\ now solution of the amplitude 
equation (QlJ The dimcnsionlcss constant is a w 13.5. 
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Doping 


r,[K] 


Vq[K] 


To[K] 




0.92 


92.9 


261.9 


124.3 


13.7 


0.87 


93.7 


178.0 


143.3 


14.3 


0.80 


88.8 


121.8 


158.9 


10.6 


0.73 


69.3 


111.0 


177.6 


13.7 


0.67 


60.8 


88.9 


193.0 


18.1 


0.57 


55.9 


76.1 


213.5 


20.6 


0.48 


47.7 


57.6 


241.5 


16.0 


0.43 


31.5 


37.3 


251.2 


9.0 



TABLE I: Extracted parameters Vo,To,a from the fits of 
YBCO specific heat, (see Figure 0. 

In Fig. |S1 the experimental speeifie heat of 
YBa2Cu306.73 Q is fitted using the variational method 
reproducing the double peak structure: a sharp peak 
below coming from phase fluctuations and a wide 
hump below T* rounded by amplitude fluctuations. The 
crossover temperature T^, where phases become random, 
corresponds to the temperature where 71 is less than 
approximately 2% of the normal specific heat. In am- 
plitude equation ijl}, a 2 dimensional density of states 
D{e) = l/W is used with W ^ 5000K, ^J. = 0.25M^ 
and U = 959K. These parameters gives Tq k, 200K and 
-00 ~ 2.14To in agreement with experiments ^3|- The 
other parameters are Vq = 72K and 77 = 5. 



B. Specific heat for the average value method 
(d-wave) 

Following equation (|19(l . the specific heat is given by 

7 = 7o((|V'I>gl)+7i (26) 

where 71 is divided by Tc instead of T since Sqi^ is a 
classical action and does not satisfy the Nernst theorem. 
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FIG. 7: YBa2Cu306+i: specific heat for different oxygen dop- 
ings X compared to the average value method. 
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FIG. 8: Average amplitude corresponding to YBa2Cu306+a: 
specific heat of Fig. Qfor difi'erent oxygen dopings x. The am- 
plitude for doping x ~ 0.9 resembles more to BCS behaviour 
whereas the amplitude for underdoped system is always very 
large due to amplitude fluctuations and large value of Tq. The 
pseudogap energy scale Eg is shown for doping x = 0.48. 
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FIG. 9: Reduced average amplitude corresponding to 
YBa2Cu306+a; specific heat for different oxygen dopings x. 
The dashed lines are average amplitudes calculated for a com- 
pletely disordered phase field. Tphi is the temperature where 
the average disordered amplitude deviates from the average 
amplitude, i.e. when effects from phases are visible. Dopings 
are the same as in figure |H| 



The average value method is compared to specific heat 
obtained for different doping in Fig. 13 For underdoped 
systems x < 0.80 we use simulations in d = 2. For the 
more overdoped, x > 0.80, simulations are done in = 3. 
The parameter 77 is fixed to 3. Parameters Vq and Tq ex- 
tracted from the fits are shown in the phase diagram of 
Fig. El and in table ^ Values of the phase specific heat 
normalisation constant a range from 10 to 20. The fitting 
procedure is done completely automatically by a random 
walk in the parameter space {Vo,To,a}, until the error 
between experimental data and the fit is minimal. As 
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FIG. 10: Underdoped La2-a;Sra;Cu04 specific heat for dif- 
ferent strontium dopings x compared to the d-wave average 
value method. 

usual, a local minimum can be reached by using this pro- 
cedure, and one can be trapped in this minimum. How- 
ever, when the dimension of the parameter space is 3 as 
in our case, different local minima can be easily excluded 
and the best fit can be achieved. 

The temperature is derived from figure |^ where the 
average amplitude is compared to the average amplitude 
(IV'Dr computed for a field where phases are random, i.e. 
the gradient of phases is set to 7r/2: 

(27) 

where D is the dimension of the system and the normal- 
isation factor Zr is 

Zr = j 1^1 e-*['''(*l'^l'+il^l')+^|-^'l'] (28) 



C. Specific heat for s-wave symmetry 

The measured specific heat of YBa2Cu306.76 Q is 
compared to our results for the s-wave symmetry in Fig. 
Illlbv using the average value method. 
We took the following values for the two parameters: 
To = 2-ibK is of the order of T*, Vb = lOSif is of the 
order of Tc- For the size of we took the BCS d-wave 
relation: il){T = 0) = 2.14To since it seems to be more in 
agreement with experiments |l7j than the BCS relation 
of i]:{T = 0) = 1.76To. This choice plays only a little 
quantitative role. 

The mean amplitude, standard deviation and phase 
stiffness are shown in Fig. E| for the same parameters 
as in Fig. The coherence length ^ is of the order the 
lattice constant when the temperature is reached. 

Considering Fig. 1111 we see that an s-wave computa- 
tion fits equally well experiments as the d-wave in Fig. 
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FIG. 11: Measurements (points) of the specific heat of 
YBa2Cu306.76 divided by "ynT. The total s-wave specific heat 
(thick blue) is the sum of the critical XY contribution (thin 
blue) and the amplitude contribution (dashed blue). If the 
amplitude of the gap were constant of size 2.14To, the contri- 
bution would be the dotted black line. 
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FIG. 12: The average amplitude {\'4>\) (thick) is large even 
above Tc whereas the coherence length ^ vanishes rapidely 
above Tc. The standard deviation - (dashed 

line) is large and almost constant. The phase stiffness in x 
direction (dotted red) jumps from 2/n to zero at Tc. Lattice 
size: iV = 80^ Parameters are To = 235 K, Vo = 108 K, 
rf = 8. 



Hence, we cannot decide which symmetry is favored 
by looking at specific heat data, s-wavc and c?-wave sym- 
metries are essentially different at low temperatures, d- 
wave specific heat has an algebraic increase with tem- 
perature whereas s-wave specific heat has an exponential 
behaviour. Therefore the s-wave contribution is smaller 
at low temperature than the d-wave. However, a lower 
s-wave contribution can be compensated by changing pa- 
rameters To and Vq. 



8 



DIFFERENTIAL CONDUCTANCE 



VI. EXTRACTED PHASE DIAGRAM 



The differential conductance between a normal metal 
and a superconductor is directly related to the density 
of states and the amplitude of the pairing field by the 
standard formula: 



(29) 



N{0) d{eV) 

where Gnn is the differential conductance between two 
normal metals. The reduced energy is ^ = e — fi. The 
s-wave density of states is according to BCS theory: 



Ofcourse, iV,(O=0if|C|<(|V^|). 



(30) 
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FIG. 13: Tlie differential conductance of equation (1291 shows 
a pseudogap, i.e a partial suppression of the density of states, 
up to temperatures near T* « 300. Inset: the differential con- 
ductance (green curve) atV = recovers its normal behaviour 
near T* . The low temperature approximation of equation II31II 
is shown by the dashed line. 

In Fig. 1131 the s-wave differential conductance nor- 
malised with Gnn is presented using an average ampli- 
tude with the same parameters as in Figure 4 of refer- 
ence (To = 159.2 K, Vo = 59.3) for temperatures from 
20 K to 300 K with 40 K intervals. One can see that 
the width of the differential conductance is proportion- 
nal to the amplitude {\ip\); and that the gap fills up due to 
thermal energy. Indeed, the low temperature differential 
conductance at y = is approximately given by: 



dl 



(1/ = 0) « Gn 



2^(1^1) ,-{|V'l)/T 



(31) 



Theses results are in agreement with scanning tunnel- 
ing microscopy on Bi2Sr2CaCu208 of Rcnner et al 
where it is observed that the pseudogap gradually fills up 
whereas its width remains constant. 

Do phase fluctuations contribute to the pseudogap? 
Yes but only for T < where the coherence length ^ 
is of the order of the lattice spacing as seen in Fig. ^1 
As shown by Eckl et al , a constant amplitude of size 
2To maintains a pseudogap in the density of states up to 
T* whereas phase effects disappear near Tc. 



Values for Tq, Vq and extracted for the specific heat 
in Fig. [Tlare reported in Fig. 1141 is obtained by using 
the comparison between amplitude from simulations and 
amplitude in the disordered phase as shown in Fig. |3 
Tff, is defined as the temperature where (|V'l)/(lV'l)r = 
h for various thresholds h = 2%,3%, .... The random 
amplitude {\ilj\)r comes from equation (|?7jl . 
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FIG. 14: Extracted phase diagram of YBa2Cu306+a:. Effects 
of amplitude fluctuations are large in the quasi uncorrelated 
pseudogap region (below To) of the copper oxides phase dia- 
gram, whereas phase correlations remain important only be- 
low r^4. The temperature T* , where observables cross over to 
normal behaviour, is located in the hatched area. 



It is remarkable that phase correlations above grow 
rapidely in the underdoped regime following the line, 
and reduce when approaching the overdoped regime. 
lines are similar to Nernst effect results @. is also 
in agreement with the evolution of the temperature Tscf 
(scf=superconducting fluctuations) of S. H. Naqib et al 
|2| derived by transport measurements. The gradient 
specific heat seems to disappears more rapidely. This 
latter doping dependence is in better agreement with 
the phase diagram derived in Hall effect experiments |20l | . 

In fig. 1141 the pseudogap energy scale Eg is reported 
in the phase diagram of YBCO. The pseudogap energy 
scale Eg is defined here as the amplitude at T = 200K: 



E„ 



m)T-. 



:200K 



(32) 



Eg shows approximately the same doping dependence as 



the one found by Loram and Tallon |J|. However Eg is 
not related to some additional hidden order parameter 
since Eg is simply related to amplitude fiuctuations as 
shown by calculations. 
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FIG. 15: The pseudogap energy scale Eg in the extracted 
phase diagram of YBa2Cu306+a;. 

VII. DISCUSSION 

We have show that amphtude and phase fluctuations 
can explain the emergence of the pseudogap region of 
underdoped high temperature superconductors. Phase 
coherence disappears completely near a temperature 
above Tc, and therefore, for T > T^, the pseudogap 
region is dominated by amplitude fluctuations. We 
flnd that the mean field temperature Tq has a similar 
doping dependence as T*, signaling that the pseudogap 
region is due to independent fluctuating pairs. The 
energy scale Eg of the pseudogap is also derived. The 
large value of Eg in the underdoped domain is due to a 



finite To and large amplitude fluctuations. At maximum 
doping, To remains of the order of Tc but amplitude 
fluctuations are much weaker. Therefore Eg vanishes 
linearly when approaching maximum doping. The 
fact that Eg takes small values above maximum doping 
should not be interpreted as a quantum critical point but 
as a crossover: amplitude fluctuations are still present 
in the overdoped regime and Eg has always a finite value. 

Comparison with measured specific heat on under- 
doped LSCO reproduces the double peak structure like 
in YBCO: a sharp peak below T^ coming from phase 
fiuctuations and a separate wide hump below T* rounded 
by the amplitude. The spin susceptibility, related to 
the amplitude, recovers its normal behaviour near T* 
whereas the orbital magnetic susceptibility, related to 
phases, disappears near T^. These considerations are 
independent of the underlying pairing mechanism, and 
any microscopic theory inducing pairing should lead to 
similar conclusions. 

An important difference between overdoped and un- 
derdoped superconductors is the "separation" between 
effects of amplitudes and phases: in the overdoped 
regime, contributions of amplitude and phase are 
added producing only one peak in the specific heat for 
example. In the underdoped regime, phase correlations 
remain near Tc and still produce a small peak whereas 
amplitude fluctuations extend to much larger tempera- 
ture producing a separate wide hump between Tc and T* . 

Further work is needed in order to include in the theory 
effects of the magnetic field on the pseudogap. 

This work has been supported by the Swiss National 
Science Foundation. 
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